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Yielding transitions in athermal amorphous materials undergoing steady-state shear flow resemble 
critical phenomena. Historically, they have been described by the Herschel-Bulkley rheological 
formula, which implies singular behaviors at yield points. In this paper, I examine this class of 
phenomena using an elementary version of the thermodynamic shear-transformation-zone (STZ) 
theory, focusing on the role of the effective disorder temperature, and paying special attention to 
scaling and dimensional arguments. I find a wide variety of Herschel-Bulkley-like rheologies but, 
for fundamental reasons not specific to the STZ theory, conclude that the yielding transition is not 
truly critical. In particular, for realistic many-body models with short-range interactions, there is a 
correlation length that grows rapidly, but ultimately saturates near the yield point. 


I. INTRODUCTION 

Numerical simulations and analytic approximations 
imply that yielding transitions in athermal amorphous 
materials, undergoing steady-state shear flow, resemble 
critical phenomena. These transitions are characterized 
by fluctuating regions of correlated, irreversible, molec¬ 
ular or granular rearrangements, whose sizes grow as 
the shear rates decrease and the systems approach their 
yield points. Recent work on this subject has been in¬ 
fluenced by the stochastic model introduced in 1998 by 
Hebraud and Lequeux.[l:] For example, Bocquet, Colin 
and Adjari Q, and Nicolas, Martens and Barrat [3j 
have studied related models that focus on elastic inter¬ 
actions between localized, stress-driven events. Lemaitre 
and Caroli Q have studied yielding in zero-temperature, 
molecular-dynamics simulations of two-dimensional, bi¬ 
nary, Lennard-Jones mixtures. Lin et al. (Sj and Muller 
and Wyart 0 have developed new techniques for study¬ 
ing instabilties and scaling behaviors for models of non- 
crystalline materials near yield points. Common out¬ 
comes of these investigations are versions of the well 
known Herschel-Bulkley (HB) rheology in which the in¬ 
cremental stress above the yield point is proportional, 
in many of these cases, to the square root of the plas¬ 
tic shear rate. This singular behavior, if accurate, would 
imply that some kind of collective motion is occurring at 
the transition. 

From its inception, the theoretical picture that has 
come to be known as “Soft Glassy Rheology” (SGR) 0[|| 
focused explicitly on Herschel-Bulkley and related rheo¬ 
logical behaviors. In contrast, the shear-transformation- 
zone (STZ) theory 0 G3 was never presented clearly 
in rheological language, which now seems to me to have 
been unfortunate. The STZ theory is an attempt to de¬ 
scribe amorphous plasticity in terms that are closer than 
SGR to the underlying physics and statistical thermo¬ 
dynamics of these processes. It is an intrinsically mean- 
field theory; it does not explicitly describe correlations 
between localized shear transformations, at least not in 
its present form. However, it does describe mean-field 
collective motion via a thermodynamically defined, effec¬ 
tive disorder temperature. The noise generated by driven 


plastic deformation disorders the system, raises the effec¬ 
tive temperature, and thereby creates new flow defects 
(STZ’s) that, in turn, contribute to the deformation rate. 
This nonlinear mechanism provides a framework in which 
to describe a variety of Herschel-Bulkley-like rheological 
behaviors. A similar mechanism has been shown to ac¬ 
count for the basic features of dislocation-induced plastic¬ 
ity in polycrystalline materials , where the generation 
of defects, i.e. dislocations, is a hardening rather than a 
softening mechanism. 

In fact, HB behavior has appeared implicitly in several 
earlier STZ-related publications. For example, in [l2j 
(LM), Manning and I analyzed data from glass-dynamics 
simulations by Haxton and Liu [13] (HL), who saw HB 
power laws at low temperatures and in a range of rel¬ 
atively high shear rates. Importantly for present pur¬ 
poses, HL also determined effective temperatures directly 
by measuring pressure fluctuations. Similar rheological 
results appear (less clearly) in my papers with Egami 
M and Lieou [jlij]. In all three of these investigations, 
our primary interest was in the transitions from yield¬ 
ing to viscous behavior as functions of temperature or 
packing fraction. We presented our results in the form of 
log-log plots, and thus did not pay attention to the rheo¬ 
logical significance of the low-temperature data, nor did 
we pay close enough attention to the yielding transition 
itself. My purpose here is to look harder at the latter as¬ 
pects of the theory. I particularly want to emphasize the 
role of fundamental, dimensional and scaling arguments. 
It is these general arguments, and not system-specific 
ones, that lead to the conclusion that yielding in ather¬ 
mal amorphous materials is not a true critical transition. 


II. STZ BASICS 

A central feature of the STZ theory is that it treats 
the density and orientations of localized flow defects 
as dynamical variables. The STZ’s fluctuate into and 
out of existence in the environment of an elastic solid. 
When present, they undergo irreversible rearrangements 
and thus produce plastic deformation in response to 
stresses. These local transitions have been observed di- 
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rectly in simulations and, recently, in a detailed experi¬ 
mental study of a colloidal glass by Jensen, Weitz, and 
Spaepen. jl||. The STZ equations of motion that deter¬ 
mine these behaviors are subject to constraints imposed 
by the first and second laws of thermodynamics. This 
theory has been discussed extensively in the literature 
[lOl Il7| along with its applications to shear banding 18], 
fracture toughness (a, oscillatory viscoelasticity |20| |. 
and the like. In what follows, I briefly summarize the 
main features of the STZ theory without repeating de¬ 
tailed derivations. As part of this summary, however, I 
raise several fundamental issues that I think need special 
emphasis. 

For athermal systems undergoing steady-state flow in 
pure shear, the basic STZ relation between deviatoric 
stress a and plastic shear rate j pl has the form 

q = 'y pl tq = eg e~ 1/x f(a). (2.1) 

Here, \ is the effective temperature expressed in units of a 
characteristic STZ formation energy; and the Boltzmann 
factor e -1 / x is proportional to the density of STZ’s. eo is 
a dimensionless constant, roughly of the order of unity. 
It is the product of the prefactor of the Boltzmann ex¬ 
ponential, i.e. a site density, and the volume of the de¬ 
formable core of an STZ. /(cr)/r 0 is the stress-dependent 
STZ transition rate, to be specified below. For simplic¬ 
ity, I consider only pure shear in the x, y plane, so that 

q — Qxx — Qyy aa d ® — &xx — U "yy 

Note that I have introduced a time constant To in 
Eg. (12. II) . Even in athermal systems such as those con¬ 
sidered here, the internal dynamics must determine a 
nonzero time scale on which local configurational fluc¬ 
tuations relax and dissipate energy. By introducing To, 
I emphasize that I am not considering athermal, quasi¬ 
stationary, numerical models, which exclude any notion 
of time whatsoever, and therefore do not allow us to dis¬ 
tinguish between fast and slow processes as is done im¬ 
plicitly in the definition of q in Eq. m■ Nor am I con¬ 
sidering models with infinitely long ranged interactions, 
or even jammed granular materials in which force chains 
may span the system, tq must be a well defined, few-body 
relaxation time. Systems that are driven faster than 1/ro 
behave differently than those that are driven more slowly; 
thus q ~ 1 plays a special role in the following analysis. 

Derivations of the Boltzmann factor on the rig ht-hand 
side of Eq. EU appear in STZ papers such as PL), H3] • As 
shown there, this elementary statistical formula is valid, 
not just for equilibrated systems, but also for systems 
that are driven persistently away from equilibrium. It fol¬ 
lows directly from the requirement that the statistically 
defined entropy be a non-decreasing function of time. \ is 
the derivative of the configurational energy with respect 
to the configurational entropy (see Sec. EU; it should not 
be thought of as a “noise temperature” or used, for ex¬ 
ample, in a Langevin equation - certainly not without 
some systematic rationale. 

For present purposes, the rate factor in Eq. EU can 


be written in the form 

/(a) = ( C ^( 1 - Ft) if H> 

1 0 if |oj < 

The yield stress a y locates an exchange 
bility between non-flowing and flowing states. The fully 
time dependent STZ equations of motion (see 0 ) de¬ 
scribe the behavior of the density of STZ’s and their 
average orientations with respect to the stress. These 
equations have stable fixed points for flowing states only 
if |oj > <T y , where the rate at which STZ’s are deacti¬ 
vated by making transitions into the forward direction is 
balanced by the rate at which new STZ’s are formed. In 
the athermal limit, cr y /\cr\ is the fraction of STZ’s that 
is aligned with the stress; in Eq. (12.21) . it plays the role 
of a back stress. a y is a system-specific quantity. It first 
appears as a factor in the relation between the rate of 
work done by the external driving force and the strength 
of the mechanical noise that induces STZ creation and 
annihilation. (See the discussion of the noise frequancy 
T in Sec |V]) 

The factor C(tr) in Eg. (12.2D is a linear combination of 
STZ transition rates, and is a symmetric function of a. 
In LM, Manning and I wrote this term in the form 


Uyi ( 2 . 2 ) 
a y . 

of dynamic sta- 


C(cr) 



(2.3) 


Here, I assume that any overall multiplicative constant 
on the right-hand side of this equation has been incorpo¬ 
rated into the rate factor l/r 0 . For larger stresses, of the 
order of or greater than some ay, Manning and I assumed 
that this rate would grow as the n’th power of the stress. 

There are numerous possible rationales for choosing n. 
For example, n = 1 would be appropriate for a colloidal 
suspension in which particle motions are subject to linear 
viscous drag. A different possibility, for systems in which 
particles interact via short-ranged repulsive forces, is to 
invoke Bagnold scaling. As conjectured in LM, if there 
are no natural stress scales in the system other than the 
rate-independent cr y and ay, then dimensional analysis 
might require that the dynamic stress a be proportional 
to an acceleration, i.e. an inverse time squared, and thus 
be proportional to the square of a dynamic rate. Con¬ 
versely, the rate C(ct)/to would then be proportional to 
the square root of the stress. This result, by itself, is in¬ 
consistent with HB rheology. Nevertheless, the Bagnold 
assumption with n = 1/2 fits the Haxton-Liu data quite 
nicely because, as shown in LM, the effective tempera¬ 
ture x and the corresponding STZ density are functions 
of q = 7 p/ ro, which brings the time scale To back into 
the theory. This Bagnold analysis seems to me to be 
plausible but not compelling in the present context. 
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III. EFFECTIVE DISORDER TEMPERATURE 

The equation of motion for the effective temperature 
X is a statement of the first law of thermodynamics, i.e. 
energy conservation. Let U denote the configurational 
internal energy as a function of the configurational en¬ 
tropy S. Then x = dU/dS ; and the first law takes the 
form 


to \ S « r 0 V c e ff X = 2 V q a - Q, (3.1) 


where x S is the rate of change of U, c e ff is an effective 
specific heat, 2qa/ro is the rate at which work is done 
on the system per unit volume, V is the volume, and 
Q/tq is the rate at which heat is dissipated. The second 
law requires that Q = V k(\ — 0), where k is a non¬ 
negative heat transfer coefficient, and 9 is the ambient 
temperature expressed in the same units as those used 
for x • For the athermal systems considered here, 6 sa 0 . 
Thus, the steady-state equation of motion for x can be 
written in the form 


toX 


2qa 

C eff 


~ «X = 


2qa 

Ceff 


1 - 


X 


x(q). 


= o. 


(3.2) 


With these assumptions, I have defined x — x(q) bi 
steady state. Note that the underlying physics is thermal 
transport, which somehow must be encoded in x(q). 

To evaluate x(q), it is useful to think about the anal¬ 
ogous relation between ordinary temperature and relax¬ 
ation rates in glasses. Note first that, if we shear or 
otherwise “stir” a system more slowly than any of its in¬ 
ternal relaxation rates, i.e. set 5 < 1 , then this system 
must ultimately reach a g-independent effective temper¬ 
ature, say xo- Shearing an amorphous material necessar¬ 
ily involves rearranging the particles, i.e. forcing them 
to move around each other to accommodate the shear 
deformation. As we shear slowly, we bring clusters of 
particles into unstable configurations, from which they 
relax relatively quickly into new stable positions. These 
persistent rearrangements produce a fluctuating steady 
state of disorder, i.e. an effective temperature Xi that is 
independent of the shear rate so long as that rate is slow 
enough. Achieving the steady-state value of x requires 
only that the accumulated shear strain be large enough, 
perhaps of the order of unity, so that a statistically signif¬ 
icant number of rearrangements occurs. This argument 
should be valid whether or not the STZ model of those 
rearrangements is accurate. It is clearly consistent with 
the HL molecular dynamics simulations shown in Fig. 2. 

This rationale for asserting that x _• v Xo > 0 in the 
limit q —> 0 is crucial for understanding the fundamen¬ 
tal distinction between realistic many-body systems and 
some of the models that have been proposed to describe 
yielding, e.g. refs. [1-3] or the recent AQS study by 
Salerno and Robbins plj . The authors of [21] cannot 
compare the imposed strain rate with dynamic relax¬ 
ation rates, as is necessary in my argument about xo- 
They also have no way to determine whether extended 


cascades of events might be affected - perhaps limited - 
by other events that can occur simultaneously. Those 
questions are also well beyond the scope of the simpler 
models described in refs. [lJ-Q. This fundamental obser¬ 
vation leads me to conclude that, for purposes of study¬ 
ing yielding transitions, these models are not in the same 
universality class as realistic many-body systems. 

The glassy analog of this behavior is that the configu¬ 
rational relaxation time r a , like 1 /q, diverges (or becomes 
unmeasurably large) as the system approaches an ideal 
glass transition. Conversely, the temperature varies in¬ 
creasingly slowly as r a goes to infinity. Like the glass 
transition temperature, xo is a system-specific parame¬ 
ter. We will see in Sec. m that its nonzero value deter¬ 
mines the nature of the yielding transition. 

At the other end of the range of shear rates, the 
glass analogy implies Arrhenius behavior. Sufficiently 
far above the glass temperature, the configurational re¬ 
laxation rate t~ x (or the inverse viscosity, or the dif¬ 
fusion constant, or the like) appears experimentally to 
be controlled by a thermally activated process with a 
temperature-independent energy barrier. Similarly, for 
values of q not too small, the Haxton-Liu measurements 
of the effective temperature are well described by a rela¬ 
tion of the form q ss q 0 e~ A ^ x , where A is an activation 
energy expressed here, like x> in units of the STZ forma¬ 
tion energy. To interpolate between these two limiting 
behaviors in [l^, Uj, > m Y coauthors and I have postu¬ 
lated a modified Vogel-Fulcher-Tamann (VFT) formula: 


q(x) = 9o exp 


— ~-Oteff(X ) , 
X, 


(3.3) 


where 


a eff(x) 


Xi 


X ~ Xo 


exp 


b X ~ Xo ' 
XA - Xo. 


(3.4) 


The first term on the right-hand side of Eg. (13.411 is the 
conventional VFT divergence at x = Xo- Here it is cut 
off in x by an exponential factor so that the behavior of 
q(x) in Eg. (13.31) is dominated by the Arrhenius term for 
X — Xo > Xa — Xo- This interpolation formula has no 
special physical significance that I am aware of and, so 
far, its details have not seemed crucial for interpreting 
experimental or computational data. On the other hand, 
the parameter A plays a special role in what follows. 

In interpreting Ea. (13.3D . note that the STZ theory has 
an absolute upper limit of validity where x —>■ 00 . At such 
large values of x, the density of STZ’s becomes large, and 
the theory is no longer consistent with a model of a solid 
containing a dilute population of flow defects. A natural 
way to formulate the theory, then, is to choose q 0 = 1 so 
that the maximum theoretical shear rate is = 1/to- 
Above this rate, the system has insufficient time to relax 
between irreversible rearrangements; i.e., it behaves like 
a fluid. Adopting this convention, however, means that 
we need to use a physically realistic estimate of tq . 
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FIG. 1: Log-Log plot of stress a in units of the yield stress a y 
as a function of dimensionless shear rate q. The data points 
are from HL, with rescaled values of q as described in the text. 



FIG. 2: Semi-log plot of the steady-state effective tempera¬ 
ture x as a function of dimensionless shear rate q. The data 
points are from HL, with rescaled values of q as described in 
the text. 


IV. HERSCHEL-BULKLEY BEHAVIOR 


At large stresses and shear rates, the preceding set of 
equations immediately produces a generalized Herscliel- 
Bulkley relation. In this limit, 


<7 


e 0 e 


-i/x 



Qoe 


-A/X' 


Eliminating y, we find 



P = 


A- 1 
nA 


(4.1) 


(4.2) 


Then, for example, the Bagnold choice n = 1/2, and A = 
4/3, produces the conventional HB power law, j3 = 1/2. 
However, with different choices of parameters, this STZ- 
based theory of athermal yielding produces a wide range 
of rheological behaviors, well beyond the simple square- 
root law. 

To explore some of these possibilities, I start by re¬ 
analyzing the low temperature HL data. I have made 
two changes from the earlier LM analysis. First, instead 
of using the microscopic time scale adopted by HL, I 
have deduced a presumably more realistic to by setting 
go = 1- The resulting r 0 is about a factor 12 larger than 
the original HL value, but is still a microscopic time scale. 
Second, in our previous attempt to fit all the HL data, 
including at high temperatures, with as few parameters 
as possible, Manning and I used what I now think was 
an unrealistically small value of ay, and thus needed to 


use a correspondingly small value of eo. In simple ather¬ 
mal systems, however, there is really only one physically 
meaningful stress scale, which is set by the shear modulus 
/r. Both cr y and ay should be very roughly of the order 
of /r. Therefore, I have set ay = a y , and have measured 
the dynamic stress a in units of a y . Then I have ad¬ 
justed eo to fit the data, finding eo = 0.26 in accord with 
my expectation that this parameter should be roughly 
equal to unity. Other parameters, as given in LM, are: 
A = 1.5 (so that j3 = 2/3 as observed), \o = 0.2 (directly 
observed as shown in Fig. |2|), Xl = 0.26, XT = 0.3, and 
6 = 3. 

Figures [I] and [ 2 ] show the HL data in their original 
form as a log-log plot of cr/a y versus shear rate g, and a 
semilog plot of the effective temperature x(g), along with 
the theoretical curves evaluated using the full equations 
given above. In Fig. [2j note that the upper limit of the 
data at q = 0.25 occurs at x — 0.7, which should be about 
at the limit of validity of the STZ theory, as intended. 
Figure[3]is a direct plot of the data in which the HB form 
is easily visible, with the crossover to the q <C 1 behavior 
squeezed into a small part of the graph near q = 0. For 
comparison, the dashed curve shows the approximation 
in Eo. (14.21) displaced upward by one unit of the yield 
stress. As can just be seen in both the data and the 
theory, the stress rises linearly above the yield point: 


a (Jy 


q — e 1 /* 0 . 

eo 


(4.3) 


For the chosen parameters, this formula means that the 
initial slope of the graph in Fig. [3] is large but finite, 
approximately 570. The rapid growth of x and the cor- 
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FIG. 3: Stress u in units of the yield stress a y as a function 
of dimensionless shear rate q. The data points are from HL, 
with rescaled values of q as described in the text. The dashed 
line is the approximation given in Ea. (l4.2l) displaced upward 
by one unit of the yield stress. 


FIG. 4: Stress u in units of the yield stress a y as a function 
of dimensionless shear rate q, for five different values of the 
dimensionless activation energy A. From bottom to top, A = 
2.0, 1.33, 1.1, 1.0, and 0.9. The corresponding values of the 
HB exponent are (3 = 1.0, 0.5, 0.18, 0, and — 0.22. 


responding growth of the STZ density cause the material 
to soften, with the result that the curve bends over into 
power-law behavior at small values of q. 

To illustrate the variety of rheological behaviors that 
emerge from this theory, I show in Fig. |4] a set of five 
different curves of stress versus shear rate, using all but 
one of the same parameters that were used to fit the 
HL data shown in Figs. |T) - [3] The exception is that 
I have chosen a sequence of different values of the di¬ 
mensionless activation energy A that controls the dissi¬ 
pation rate in Eg. (13. 21) . From bottom to top in the fig¬ 
ure, these values are A = 2.0, 1.33, 1.1, 1.0, and 0.9. 
The corresponding values of the HB exponent are (3 = 
1.0, 0.5, 0.18, 0, and — 0.22. The second curve, with 
/3 = 0.5, is the square-root law that appears in the 
fluctuation-theory literature, e.g. 0,1 a ; and the third 
and fourth curves, with small values of /3, are included 
just to show the transition between stable and unstable 
rheologies. 

The more interesting cases shown in Fig.[4]are the first 
and the fifth. The first, at the bottom with /3 = 1, illus¬ 
trates one of the ways in which this theory can produce a 
conventional, linear, Bingham rheology. By far the most 
realistic way for the linear law to be observed, however, 
is when tq is a microscopic time, perhaps a few molecu¬ 
lar vibration periods, so that shear rates with q ~ 1 are 
out of the range of most experiments. Then, Eg. (14.3D . 
with a constant xo, is accurate throughout an observable 
range of shear rates, q -C 1. This is what happened, 
for example, in earlier analyses of plasticity in metallic 
glasses (e.g. HU), where Xo was more nearly of the order 


of 0.1 instead of 0.2. Yet another possibility, still within 
the scope of this theory, is when materials are intrinsi¬ 
cally disordered and soft, i.e., when xo is large, so that 
the initial slope of cr(q) in Eg. (14.3D is small and remains 
linear at observably large shear rates. 

The fifth interesting case in Fig. [4] at the top with 
/3 = —0.22, exhibits shear-rate weakening above q = 0.07, 
and therefore must be dynamically unstable. Daub and 
Carlson [23[ have examined a model of this kind in detail, 
and have shown that it produces both shear-banding and 
stick-slip instabilies. More recently, Lieou et al. [24j have 
applied this kind of theory - but with (3 > 0 - to granu¬ 
lar materials of the kind that occur in fault gouge. They 
have added terms to k in Eo. (l3.2D to account for shear- 
rate dependent frictional dissipation due to collisions be- 
ween angular grains, and also to account for dissipation 
induced by tapping or other vibrational perturbations. 
In this generalization of the theory, the effective temper¬ 
ature can become a non-monotonic function of q , and 
the volume may also vary non-monotonically. But these 
modifications are beyond the scope of the present discus¬ 
sion, because they implicitly introduce new intrinsic time 
scales analogous to r 0 . 


V. SPATIAL HETEROGENEITIES 

As a first step in studying spatial variations in this 
theory, I have linearized the equation of motion for y, 
Eq. ( 13.2D . in the space and time dependent variable y' = 
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FIG. 5: (Color online) Correlation length £/a as a function of 
dimensionless shear rate q. The two sets of experimental data 
(blue triangles and red squares) are for the two larger forcings 
shown in Fig.5 of Jop et al [2y], and the theoretical scales of 
length a and shear rate q are chosen here to be consistent 
with that data. 


X ~ x(s)> and have added a diffusion term: 


To x' ~ 


2 a_q x/_ 

c eff x(q) 


+ VV 2 


X' 

x(q)‘ 


(5.1) 


Because x(q) diverges when q —> qo, it is best to linearize 
in x'/ x(q) i and to define the diffusion coefficient V ac¬ 
cordingly. Diffusion terms of this kind hav e ap peared 
previously in the STZ literature, e.g. in [il, 23]; but it 
seems to me that T> has never been given the attention 
that it needs. V must have the form T £ 2 , where T/r o is 
the noise-generated attempt frequency for activation and 
annihilation of STZ’s, and i is an elementary diffusion 
length, which I argue must be the average spacing be¬ 
tween STZ’s. This is in contrast to my argument in p5j . 
where I considered self diffusion of particles, and set t 
equal to the size of the STZ core. Here, I am considering 
diffusion of energy or perhaps, in the spirit of Hi , stress 
fluctuations; thus, the spacing between STZ’s seems ap¬ 
propriate. 

To implement these ideas, I use primarily dimensional 
analysis. The only scalar quantity in this system with 
dimensions of rate is the power per unit volume, a q/To- 
To convert this quantity into a frequency, divide by a 
stress times a volume (an energy), say, cro vq, and multi¬ 
ply by the only relevant volume, i.e. the volume per STZ, 
voe 1 /*. The factors Vo cancel, and subsequent athermal 
STZ analysis reveals that Go is equal to the yield stress 
Cy. Thus, T = 2 q (a/a y ) e 1 ^. Similarly, for dimension 


d, t d = Vo e 1/,% . The result is 

V = 2 a 2 q — e {1+2,d)/x , (5.2) 

Gy 

where vq = a d defines a characteristic interparticle length 
scale a. 

Static solutions of Eg. (15. II) have the form e~ r ^, where 


2 _ c e ff T> _ c e ff a 2 1 + 2/d 

2a q a y P [ x(q) _ ' 

In the large-g limit where q ss qo e~ A t x , Eg. (15.3D 

c e ff a 2 fq 0 \ 2v _ 1 + 2,/d 
4 ~ Gy v q) ’ 17 _ 2H ■ 


(5.3) 
becomes 

(5.4) 


This growing length scale saturates near the yielding 
transition, where q —> 0 , and X —» xo- 

The full function £(q) is shown in Fig. [5] in comparison 
with data taken from Fig. 5 of Jop et al. [26[. The data 
shown here were obtained by measuring velocity-velocity 
correlations in a dense emulsion flowing under pressure 
through a gap between parallel glass plates. Correlations 
between velocity fluctuations perpendicular to the flow 
were measured as functions of two-point separations par¬ 
allel to the flow, at fixed distances from the center of the 
gap and, thus, at fixed shear rates. In computing the 
theoretical curve, I have used qo = 18 in order to match 
the scale of shear rates 7 reported by those authors; and I 
have set \Jc e ff /a y = 0.08 in Eq. (ED, again to match the 
units of length used in [26(. To optimize the fit, I chose 
A = 1.2, which means that the HB exponent is /3 = 1/3, 
and v = 0.69 in Ea. (l5.JD . Otherwise, the STZ parame¬ 
ters (xoj Xii XAi and b) are the same as those used in 
fitting the HL data in Figs. CD -[3] This is the best I can 
do in the absence of more detailed information about the 
rheological properties of this particular emulsion. 

I emphasize that I have not attempted here to con¬ 
struct a complete theory of the experimental observations 
in |26], even within the assumptions that led to Ea. (l5.3D . 
To do that, I would have had to couple the equation of 
motion for y, Eq. & to a position-dependent and fully 
tensorial version of the STZ flow law, Eq. (El. More 
importantly, I would have had to introduce a nonlocal 
version of that flow law in which the strain rate on the 
left-hand side is an exponentially weighted average of 
the driving term on the right-hand side within a neigh¬ 
borhood of size £. In other words, I would have had 
to develop an effective-temperature version of the nonlo¬ 
cal rheology described recently by Hennan and Kamrin 
[H [jsl] in their theory of dense granular flows. Their 
analysis must be closely related to the one presented here; 
it should be useful to explore that relationship. 

A nonlocal analysis should not make qualitative 
changes in the results shown in Fig. [5] which relate to 
measurements at fixed distances away from the center 
of the gap where shear rates are constant and nonzero. 
However, a fully nonlocal analysis would be essential for 
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computing the velocity profile across the flow, as mea¬ 
sured by Jop et al. [26|. We know by symmetry that, in 
plug flow of this kind, both the shear rate and the shear 
stress change sign from one side of the gap to the other. 
Therefore, there is a region near the center of the gap in 
which | a \ < a y , where a local relation of the kind shown 
in Eq. m cannot be valid. A nonlocal calculation of 
the velocity profile is reported in [26|, where it is used 
to estimate a length scale whose slow shear-rate depen¬ 
dence may be consistent with v = 1/4 as predicted in 
[1]. It remains to be seen whether that behavior can be 
distinguished from the saturation effect predicted here. 


VI. REMARKS AND QUESTIONS 

A principal conclusion of this analysis is that the ther¬ 
modynamic STZ theory, like SGR, predicts a wide range 
of Herschel-Bulkley-like rheological behaviors, consistent 
with a long history of experimental observations. The 
surprising result is that - largely due to the properties of 
the effective disorder temperature x ~ the yielding tran¬ 
sition in this theory is non-critical. Immediately above 
the transition, the stress rises linearly in the strain rate, 
and the correlation length remains non-diverging. This 
result contradicts almost all recent theories of yielding in 
amorphous materials. I have argued here that this dis¬ 
agreement occurs because those recent theories are based 
on models or simulation schemes that do not accurately 
describe the dynamics of realistic many-body systems. 

I emphasize that this predicted deviation from criti¬ 
cality at the yielding transition is not an artifact of any 
phenomenological assumption, or of any special feature of 
the STZ theory. On the contrary, it emerges from first- 
principles. As argued in Sec. Mil the state of disorder 
at small dimensionless shear rate q must be independent 
of q , i.e. x Xo > 0- This nonzero degree of disorder 
means that, in steady-state deformation, there must be a 
constant density of some kind of flow defects and, there¬ 
fore, a linear relation between q and the stress increment 
just above the yield point. This theoretical conclusion is 
supported by a close look at the numerical simulations 
shown in Fig. 0 and by the experimental data shown in 
Fig. [5] 

On the other hand, I emphasize that this unexpected 
result does not in any way invalidate the evidence that 
yielding transitions are generally accompanied by large 
stress fluctuations and cascades of extended, correlated 
events. There are two complementary kinds of questions 
that need to be asked in this regard. First, the exist¬ 
ing theories have been extremely useful in identifying 
the physical mechanisms that underly yielding dynam¬ 
ics. In what ways are they accurate? Where, precisely, 
might they be missing essential ingredients? Second, the 
mean-field STZ theory provides only an approximate de¬ 
scription of the interactions between spatially separated 
yielding events. In what ways might it be modified in or¬ 
der to be more accurate? Does it need to be so modified? 


Theoretical evidence in favor of criticality at yielding 
transitions has come from stochastic models of the kind 
originated by Hebraux and Lequeux [l] and extended 
in [2, Q. As argued in [ 5 J and the yielding mech¬ 
anism looks as if it should be in the same universality 
class as the depinning mechanism that produces, for ex¬ 
ample, broad distributions of slipping events on earth¬ 
quake faults. Renormalization-group analyses of depin¬ 
ning models have produced scaling relations similar to 
the Herschel-Buckley law. (For example, see Fisher et. 
al. H2].) How precise are these analogies? 

The Hebraux and Lequeux model produces a Herschel- 
Bulkley exponent of exactly 1 /2 in the limit of zero strain 
rate. Even if we assume that some limiting approxima¬ 
tion is being made that eliminates the saturation effect 
(perhaps by implicitly setting yo = 0), the question re¬ 
mains: why is ft = 1/2 in this class of models? How 
might it be modified to produce other values of /3? We 
know that other HB exponents do occur, e.g. in the 
Haxton-Liu results discussed here, and in a wide variety 
of rheological situations. According to Eq. (14.21) . this ex¬ 
ponent is determined by the activation energy A, which 
controls the heat flow. Is there any analogous physics in 
[lj? Conversely, we must ask: What physics determines 
A in the STZ theory? Under what physical circumstances 
might we find unstable dynamics with A <11 

What are the corresponding strengths and limitations 
of the thermodynamic STZ theory? We know from recent 
experience that the STZ equations of motion, when gen¬ 
eralized for variations in space and time and coupled to 
the equations of motion for elastic fields, constitute a self- 
contained dynamical theory of amorphous rheology. The 
theory does a good job of describing phenomena such as 
shear banding instabilities jTsj, brittle and ductile failure 
in the neighborhoods of crack tips 0 , stick-slip behavior 
of granular materials in earthquake faults [23, [24j] , and 
the like. In analogy to the well known equations of mo¬ 
tion for fluids, it seems reasonable to expect that these 
rheological equations of motion, when solved for large 
systems coupled to external forces and boundary condi¬ 
tions, will predict chaotic behaviors with heterogeneous 
deformations and local failures on many different length 
and time scales. Will the results of such calculations be 
consistent with those based on the depinning analogy? 
Or with experimental observations? 

A related question is whether the deterministic nature 
of the rheological theory presented here might be inad¬ 
equate for describing broad distributions of event sizes. 
The scaling analysis in [29[ starts with the assumption of 
a quenched random pinning force. Hebraux and Lequeux 
couch their analysis in terms of probability distributions 
over the values of local stresses. Is that kind of analysis 
essential? Or might the observed behaviors be results of 
deterministic chaos, as was the case in our earlier slider- 
block studies of earthquake dynamics? |3j| Large-scale 
numerical solutions of the present rheological equations 
of motion might help answer such questions. 
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